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ABSTRACT: In order to develop specific genetic markers and determine the genetic diversity of Bangladeshi native cattle (Pabna, 
Red Chittagong) and exotic breeds (Sahiwal), randomly amplified polymorphic DNA (RAPD) analysis was performed using 12 primers. 
Genomic DNA was extracted from 20 cattle (local and exotic) blood samples and extracted DNA was observed by gel electrophoresis. 
Among the random primers three were matched and found to be polymorphic. Genetic relations between cattle's were determined by 
RAPD polymorphisms from a total of 66.67%. Statistical analysis of the data, estimating the genetic distances between cattle and 
sketching the cluster trees were estimated by using MEGA 5.05 software. Comparatively highest genetic distance (0.834) was found 
between RCC-82 and SL-623. The lowest genetic distance (0.031) was observed between M-1222 and M-5730. The genetic diversity of 
Red Chittagong and Sahiwal cattle was relatively higher for a prescribed breed. Adequate diversity in performance and adaptability can 
be exploited from the study results for actual improvement accruing to conservation and development of indigenous cattle resources. 
(Key Words: RAPD-PCR, Genetic Distance, Polymorphism, Cattle, Bangladesh) 



INTRODUCTION 

Cattle are the most common type of large domesticated 
ungulates of the subfamily Bovinae and belong to the genus 
Bos. Generally, cattle are raised as livestock for meat, milk 
leather, dung and draft purposes. It has been estimated that 
there were 1.3 billion cattle in the world (Bruford et al., 
2003) whereas Bangladesh had 24 million cattle, out of 
which 6 million were dairy cattle of local and crossbreds 
(Bhuiyan, 1993). The available breeds in Bangladesh are 
Red Chittagong, local Pabna, Sahiwal cross breed, Holstein 
Frisian cross breeds etc. It is speculated that the Red 
Chittagong evolved through interbreeding among the local 
cattle on the basis of natural selection. Sahiwal are a breed 
of zebu cattle, primarily used in dairy production, from the 
Sahiwal district of Pakistan, a dry region of Punjab which 
lies along the Indian-Pakistani border (Mason, 1996). The 
Bangladeshi indigenous cattle breed (Pabna) and Red 
Chittagong breed, are of particular importance because of 
their favorable meat, milk, leather and draft quality 
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although their productivity is relatively low in general. 
Dairying is considered a strong tool to develop village 
micro economy in Bangladesh (Shamsuddin et al., 2006) to 
improve rural livelihoods and to alleviate rural poverty. But, 
like other developing countries, breed substitution and 
crossbreeding programs have been practiced 
indiscriminately in Bangladesh for many years with exotic 
breeds and as a result valuable genetic resources are being 
lost or diluted. 

Genetic diversity or similarity can be measured through 
genetic markers. These have been used to determine 
evolutionary relationships within and between species, 
genera or higher taxonomic categories (Cornuet et al., 1999). 
RAPD markers, generated by the Polymerase Chain 
Reaction (PCR) have been widely used since the 1990's to 
assess intra-specific genetic variation at the nuclear level 
(Garrigan et al., 2002). This procedure detects nucleotide 
sequence polymorphisms in a DNA amplification-based 
assay using only a single primer of an arbitrary nucleotide 
sequence. The polymorphism between individuals results 
from sequence differences in one or both of the primer 
binding sites and is visible as the presence or absence of a 
particular RAPD band. Therefore, in order to determine the 
genetic variability of selective native, exotic and cross cattle 
populations of Bangladesh, the research was undertaken 
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using RAPD-PCR. 

MATERIALS AND METHODS 

Animal selection, blood sample collection and extraction 
of genomic DNA 

The experiment was carried out at the Animal Cell 
Culture and Molecular Biology Laboratory, Biotechnology 
and Genetic Engineering Discipline, Khulna University, 
Khulna-9208, Bangladesh during the period of November 
2010 to May 201 1 . A total of 20 cattle were selected for this 
study, in which five were native Pabna (M-1222, M-5730, 
M-1893, LF-1414, M-1980), ten native Red Chittagong 
(RCC-5730, RCC-178, RCC-175, RCC-609, RCC-80, 
RCC-173, RCC-81, RCC-691, RCC-82, RCC-77) and five 
exotic Sahiwal (SL-571, SL-614, SL-623, SL-637, SL-638) 
cattle from Bangladesh Livestock Research Institute, Dhaka. 
Three ml of blood was collected from each animal in a 5 ml 
EDTA vaccutainer. The blood was gently mixed with EDTA 
(present in the vacationer) and kept at -20°C until the study 
start. Genomic DNA was extracted from blood samples 
following the method described by Roe et al. (1996) with 
some modifications. About 800 |il of blood was taken in a 
microcentrifuge tube and the cell membrane disrupted with 
IX Standard Saline Citrate (SSC) buffer. Then the protein 
was precipitated by sodium acetate and digested with 
proteinase K at 55 °C for 1 h. After that, the aqueous layer 
was removed by Phenol:Chloroform:Isoamyl alcohol 
(25:24:1) and precipitated with 100% ethanol. DNA was 
reprecipitated and washed with Na acetate and 70% ice cold 
ethanol. Finally DNA samples were dried and dissolved in 
TE buffer and stored at-20°C for further use. 

Gel documentation, DNA quantification and PCR 
amplification 

Presence of extracted DNA was confirmed by gel 
electrophoresis separation on 1% agarose gel and staining 
with ethidium bromide. Ultraviolet (UV) spectrophotometer 
set at 260 nm was used to quantify the DNA concentration 
prior to PCR. The PCR amplification conditions were set 
according to Williams et al. (1990). Primers were selected 
mainly based on GC content, intensity of bands, presence of 
smearing, consistency within individuals and potential for 
population discrimination. Primers that have GC content 
more than 60% are suitable for RAPD analysis. We selected 
twelve subsets of primers supplied by Bioneer, South Korea. 
PCR was carried out in MyGenie 96 thermal block, Bioneer, 
South Korea. PCR amplification followed the method of 
Simon et al. (2007). The PCR reaction was carried out in a 
final volume of 20 |il reaction mixture containing 2 |il of 
template (4ng DNA), 2 |il of 10x reaction buffer, 2 |il of 
dNTPs (250 mM each), 2 units of Taq DNA polymerase, 1 
|il (100 picomol) of primer and 11 |il of sterile distilled 



water, respectively. The mixture was over layered with a 
drop of mineral oil. The PCR amplification conditions 
were: PCR lead at 105°C, hot start at 94°C for 4 min's, 
denaturation at 94 °C for 1 min, annealing at 35°C for 1 min, 
extension at 72°C for 2 min's, 35 cycles, final extension at 
72°C for 10 min, holding at 4°C for 10 min and end. 

Agarose gel electrophoresis and data analysis 

Amplified PCR products were electrophoresised 
(Agaro-power system, Bioneer Corporation, South Korea) 
on an agarose gel (2%) in TAE buffer and visualized by 
staining with ethidium bromide solution. Polymorphic 
bands were observed by high performance UV trans - 
illuminator (UVP-UK) and pictures were taken by GEL 
DOC/Bio-Imaging System affiliated digital canon camera. 
In the study of genetic diversity, each germplasm was 
scored manually for presence or absence of a particular 
amplification product. Data was analyzed by the software 
GelQuest and MEGA5.05 for phylogenic analysis. RAPD 
matrix of the 20 cattle samples using OPA07, OPA03 and 
OPA09 primers were assembled for statistical analysis and 
the sizes of the DNA fragments were estimated using 100 
bp ladders by co-electrophoresis with the PCR product 
DNA. The un-weighted pair group method with arithmetic 
averages (UPGMA) based dendrogram was constructed 
using the MEGA5.05. Numerical values were considered to 
construct the dendogram and presence of band (1) and 
absence of band (0) value considered for RAPD analysis. 
All possible pair-wise genetic distances were calculated 
according to Nei and Li (1979) using the following formula: 
d xy = l-(2n xy /(n x +n y )), Where, n x and n y are the numbers of 
bands amplified in individuals x and y respectively and 2n xy 
is the number of bands shared by those individuals In this 
analysis, smaller numbers are associated with more 
genetically similar individuals, whereas larger numbers 
suggest genetically dissimilarity. 

RESULT AND DISCUSSION 

The experimental cattle showed variability in genetical 
characteristics and the results indicated the effectiveness of 
RAPD analysis in detecting the level of polymorphism. 
Randomly 12 primers were tested among them three 
primers (Table 1) were matched and found polymorphic 
(Figures 1, 2 and 3). 

From the figures, it was observed that primer OPA-03 
and OPA-07 were more reproducible and produced a more 

Table 1. List of primers those showed polymorphism 
Primer no. Sequence (5'-3') (G+C) Contents % 

OPA-03 TGCCTCGCACCA 66J 
OPA-07 GGTGACGCAG 70 
OPA-09 GTTGCCACCC 70 



Khatun et al. (2012) Asian-Aust. J. Anim. Sci. 25:751-757 



753 



distinct RAPD profile than primer OPA-09. Three different 
primers generated various banding patterns and out of the 
21 bands scored, 14 bands (66.67%) were found to be 
polymorphic and 7 bands (33.33%) were found to be 
monomorphic in nature. Different RAPD patterns generated 
using the primers in the present experiment are shown in 
Figure 1 to 3. 

Table 2 shows that the frequencies of polymorphic 



bands varied from primer to primer in which primer OPA- 
03, OPA-07 and OPA-09 produced 80%, 42.85% and 75% 
polymorphic loci respectively. Primer OPA-03 produced 10 
bands where one at 400 bp and the other 700 bp were 
monomorphic in all samples. The polymorphic fragments 
for primer OPA-03 ranged from 900 bp to 2,000 bp. 
Furthermore, primer OPA-07 and OPA-09 produced 7 and 4 
bands, respectively. The range of monomorphic bands for 




Figure 1. RAPD profiles using primer OPA-03. The number corresponds to the serial of the genotypes; M is a 100 bp DNA ladder, 'a' 
and 'b' bands were monomorphic for all the experimental animals. 




Figure 2. RAPD profiles using primer OPA-07. On the gel 'a', 'b', 'c' and'd' bands were monomorphic for all experimental animals. 




2000 hp 
■ (1600 bp) 
1200 bp 
800 bp 



Figure 3. RAPD profiles using primer OPA-09. 'a' is monomorphic for all experimental animals. 



Table 2. List of primers and their monomorphic and polymorphic bands 





No. of 
bands 


No. of 


Size of 


No. of 


Size of 




Primer code 


monomorphic 


monomorphic 


polymorphic 


polymorphic bands 


Polymorphism (%) 




bands 


bands (bp) 


bands 


(bp) 




OPA-03 


10 


2 


700-400 


8 


2,000-900 


80 


OPA-07 


7 


4 


1,000-700 


3 


600-400 


42.85 


OPA-09 


4 


1 


1,600 


3 


1,200-800 


75 


Total = 3 


21 


7 


33.33% 


14 




66.67 
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primer OPA-07 and OPA-09 were 700 bp to 1,000 bp, 1,600 
bp, respectively. Where as, the range of polymorphic bands 
were 400 bp to 600 bp and 800 bp to 1,200 bp, respectively. 
As the number of bands common in all animals produced 
by primer OPA-07 was (4) greater than the other two 
primers (2 for OPA-03 and 1 for OPA-09), it can be 
concluded that OPA-07 may be specific for species or for 
greater zoological classes (gender, family etc.). This point 
needs to be clarified in the future by a comparison with 
other species. Sharma et al. (2009) reported, the size of 
polymorphic bands ranged between 200 bp and 2,500 bp 
using 17 operon series random primers to characterize the 
Malvi and Sahiwal breeds of cattle in Madhya Pradesh, 
India. In another study involving indigenous African cattle, 
Gwakisa et al. (1994) used RAPD markers to characterize 
the local Zebu (Bos indicus) cattle breeds of Tanzania 
(using ILO 1127 primer) and showed 61% polymorphism 
which are almost similar to the present study. The degree of 
polymorphism 66.67% is higher compared to 28% in 
taurine cattle (Talle et al., 2005) and 18% in zebu cattle 
(Marle-Koster et al., 2003). Estimation of higher level of 
genetic variation in the breeds might be consistent with the 
fact that they are highly polymorphic animals. 
Polymorphisms at the level of individual, breed, species or 
gender can be produced by RAPD-PCR but it is difficult to 
identify the bands. No breed- or sex-specific RAPD bands 
were detected in this study. Nevertheless, by mean of their 
monomorphic band profiles, some of the investigated 
primers can give the information that animals not having 
these monomorphic bands may not be related to these cattle. 
More studies should be performed on this matter in 



comparative genome studies. 

Pair wise genetic distance 

From the Table 3 it can be found that the highest value 
of Nei and Li's (1979) genetic distance (0.834) was between 
SL-623 and RCC-82. Lower estimates of genetic 
similarities indicated a high degree of genetic diversity 
between these two cattle. This may be due to their diverse 
origin and different purpose of raising these cattle. On the 
other hand, the lowest distance (0.031) existed between M- 
5730 and M-1222. It indicates, the genetic relationship 
between local cattle population was closer than others. 
Sharma et al. (2009) calculated between breed band sharing 
(BS), the value for these 2 breeds obtained was lower (0.43 
to 0.66) than within breed BS value, which was 0.70 to 0.93 
in Malvi breed and 0.68 to 0.88 in Sahiwal breed, found 
similar genetic distances of average animals like the present 
experiment. 

The present results revealed that there was a high level 
of genetic variaton among the studied cattle as indicated by 
the proportion of polymorphic loci. The values of pair-wise 
genetic distance among 20 cattle genotypes were computed 
from combined data sets for the 3 primers ranging from 
0.834 to 0.031 (Table 3). Kemp and Teale (1994) identified 
population-specific DNA polymorphisms to detect Bos 
indicus and Bos taurus in West Africa. They estimated the 
allele frequency characteristics of cattle breeds using DNA 
markers (Buitcamp et al., 1991; Mannen and Tsuji, 1993; 
Glowatzki-Mullis et al., 1995; Choy et al., 2001). 
Heterozygosity estimates obtained for the Burlina breed 
were moderate to high and consistent to the value reported 



Table 3. Estimation of pair wise genetic distance between experimental cattle populations using Nei and Li's (1979) equation 
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RCC- 


RCC- 


RCC- 


RCC- 


RCC- 


RCC- 


RCC- 


RCC- 


RCC- 


RCC- 


SL- 


SL- 


SL- 


SL- 


SL- 


M- 


LF- 


M- 


M- 


M- 


No. 


173 


178 


691 


175 


609 


77 


78 


80 


81 


82 


571 


614 


623 


637 


638 


1222 


1414 


1893 


1980 


5730 


1 

2 


0.164 
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0.111 


0.241 


*** 
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0.259 


0.172 


0.181 


*** 


































5 


0.185 


0.185 


0.393 


0.217 


































6 


0.185 


0.185 


0.304 


0.111 


0.200 
































7 
8 


0.142 
0.259 


0.142 
0.214 


0.217 
0.304 


0.111 
0.214 


0.133 
0.032 


0.696 
0.125 


0.259 




























9 


0.185 


0.214 


0.533 


0.214 


0.133 


0.307 


0.333 


0.666 


*** 
























10 


0.214 


0.142 


0.304 


0.214 


0.161 


0.096 


0.230 


0.161 


0.241 
























11 


0.071 


0.103 


0.333 


0.285 


0.803 


0.096 


0.153 


0.360 


0.172 


0.032 


*** 




















12 


0.111 


0.142 


0.185 


0.172 


0.133 


0.033 


0.142 


0.133 


0.066 


0.200 


0.333 




















13 


0.111 


0.142 


0.185 


0.714 


0.133 


0.200 


0.172 


0.133 


0.344 


0.834 


0.230 


0.200 


*** 
















14 


0.142 


0.214 


0.142 


0.142 


0.133 


0.602 


0.172 


0.062 


0.133 


0.230 


0.153 


0.307 


0.238 
















15 


0.071 


0.172 


0.143 


0.640 


0.200 


0.230 


0.241 


0.043 


0.333 


0.386 


0.214 


0.346 


0.230 


0.814 














16 


0.142 


0.103 


0.142 


0.142 


0.096 


0.386 


0.066 


0.142 


0.103 


0.161 


0.172 


0.161 


0.230 


0.083 


0.200 












17 


0.142 


0.103 


0.142 


0.103 


0.161 


0.384 


0.200 


0.142 


0.103 


0.307 


0.060 


0.307 


0.214 


0.172 


0.200 


0.360 










18 


0.185 


0.200 


0.172 


0.241 


0.161 


0.307 


0.200 


0.214 


0.096 


0.185 


0.142 


0.259 


0.172 


0.241 


0.032 


0.133 


0.033 








19 


0.103 


0.572 


0.148 


0.066 


0.217 


0.185 


0.200 


0.103 


0.062 


0.538 


0.111 


0.230 


0.083 


0.066 


0.066 


0.133 


0.142 


0.035 


*** 




20 


0.172 


0.040 


0.103 


0.133 


0.062 


0.153 


0.032 


0.160 


0.111 


0.259 


0.200 


0.259 


0.214 


0.200 


0.161 


0.031 


0.062 


0.142 


0.214 





1,2,3,4,5 



correspond to the number of the genotypes. 
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by Del Bo et al. (2001) that found an estimate of 0.68. 
However, such variability represents an important reservoir 
of genetic diversity that should be conserved. Genetic 
diversity is in fact important to allow genetic improvement 
in selected breeds as well as to facilitate rapid adaptation. 

Dendrogram analysis 

Estimation of genetic similarity within and between 
cattle and genetic distance among different cattle is an 
important application of the DNA based genetic markers. 
Dendrogram stability is an important aspect to be 
considered in genetic variability studies. An unweighted 
pair-group method of analysis (UPGMA) phylogenic tree 
based on genetic distance using the MEGA software 
(version 5.05) is shown in Figure 4. Dendogram depicted 
the relationship among twenty population of RCC, Sahiwal 
and local cattle. The dendogram segregated the twenty 
cattle into two major clusters: (I) and (II). Cluster I sub 
divided into I A and IB. Furthermore, cluster IA sub- 
clustered into IA (1) and IA (2). On the other hand, cluster 
IB subdivided into IB (1) and IB (2). 

Subcluster IA comprised of RCC-609, RCC-80, SL-637, 
SL-623, RCC-175, RCC-178, M-1222, M-5730 and RCC- 
78 cattle population and subcluster IB consisted of RCC- 
173, M-1980, RCC -691, RCC-81, M-1893, LF-1414 and 
SL-638 cattle population. On the other hand, cluster II 
belong cattle RCC-82, SL-571, SL-614 and RCC-77. From 
the dendogram, closely related indigenous population e.g., 
RCC and SL co-existed in cluster II might be due to their 



same origin as well as lower genetic variability. Whereas, 
local Pabna and RCC randomly segregated between 
subcluster IB and Cluster IA. The possible reasons for their 
random segregation could be the indiscriminate 
breeding/crossing, migration of animals, mutation in 
ancestors or due to different origin. 

Moreover, the long genetic distance between SL-623 
and local RCC-82 indicates that they can be crossed for 
parent development. The highest genetic diversity is 
alarming as it indicates that the population may deteriorate 
due to crossing with other populations of either native or 
exotic blood. In another similar study, Yu et al. (2004) used 
RAPD analysis to estimate genetic diversity and 
relationship of two native cattle breeds from the Yunnan 
province of China (De Hong cattle and Di Qing cattle) and 
four introduced beef cattle breeds (Brahman, Simmental, 
Murry Grey, and Short Horn). Using 10 primers, it was 
observed that the Yunnnan DeHong cattle breed was closely 
related to the Brahman (Bos indicus), and the Yunnan 
DiQing cattle breed was closely related to the Simmental, 
short horn, and Murry grey (Bos taurus) breeds. It is 
reported that uncontrolled hybridizations are increasing and 
the phylogenic relationship of the RCC and Local Pabna 
cattle is close because regional animal cooperative systems 
have been using RCC bulls for breeding programs since 
1995 (Ahmed, 2000). 

The high level of polymorphism obtained by the present 
study indicates that RAPD markers could be considered as 
effective tools for estimating genetic diversity in 
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Figure 4. The unweighted pair group method of analysis (UPGMA) dendrogram based on summarized data regarding differentiation 
among experimental cattle. 
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experimental cattle populations. For a more efficient 
determination of genetic diversity in cattle varieties, a 
higher number of decamer primers showing a maximum 
number of polymorphic band or other available marker 
systems could be utilized. Nevertheless, if the markers 
developed in the present experimental animals can be linked 
with quantitative trait loci for economically important traits, 
they could be used successfully in future breeding and 
germplasm conservation programs for cattle production 
systems of Bangladesh. 

CONCLUSION 

This study demonstrated that the short random 
oligonucleotides of arbitrary sequences can be used to 
successfully differentiate cattle by variation of RAPD 
fingerprints in breed-specific DNA pools and subsequently 
can be used to produce prominent fingerprints in individual 
animals. The polymorphism of the experimental cattle was 
66.67% and suggests that the gene segregation rate has 
recently increased in these animals. The genetic variation in 
the experimental cattle populations was high and therefore, 
there is scope of selective breeding (like ONBS) for future 
improvement of Red Chittagong, Sahiwal and local cattle. 
Gradual genetic erosion and extinction is threatening an 
increasing number of animal species; the major 
consequence being a loss of global genetic diversity. Sets of 
genetic markers are required to identify distinct populations 
so that they can be preserved. The local Bangladeshi Pabna 
cattle and Red Chittagong cattle are of particular 
importance because of their favourable meat, milk, leather 
and draft quality although their productivity is relatively 
low in general. Like other developing countries, breed 
substitution and crossbreeding programs with exotic breeds 
have been practiced indiscriminately in Bangladesh for 
many years and as a result valuable genetic resources are 
being lost or diluted. In addition, genetic assessment is also 
of interest for the design of genetic improvement programs 
including appropriate choice of breeds for crossbreeding. 
Thus, it is important to identify molecular characteristics of 
these breed/populations for undertaking proper breeding 
programmes and conservation strategy. 
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